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ABSTRACT 



We present seven newly discovered non-eclipsing short-period binary systems 
with low-mass companions, identified by the recently introduced BEER algo- 
rithm, applied to the publicly available 138-day photometric light curves ob- 
tained by the Kepler mission. The detection is based on the beaming effect 
(sometimes called Doppler boosting), which increases (decreases) the brightness 
of any light source approaching (receding from) the observer, enabling a predic- 
tion of the stellar Doppler radial-velocity modulation from its precise photom- 
etry. The BEER algorithm identifies the BEaming periodic modulation, with 
a combination of the well known Ellipsoidal and Reflection/heating periodic ef- 
fects, induced by short-period companions. The seven detections were confirmed 
by spectroscopic radial-velocity follow-up observations, indicating minimum sec- 
ondary masses in the range of 0.07-0.4 M Q . The discovered binaries establish for 
the first time the feasibility of the BEER algorithm as a new detection method 
for short-period non-eclipsing binaries, with the potential to detect in the near 
future non-transiting brown dwarfs secondaries, or even massive planets. 

Subject headings: methods: data analysis — planetary systems: detection — binaries: 
spectroscopic — brown dwarfs 



-3 - 



Introduction 



In a recent paper iFaigler Sz Mazeh! (120111 ) presented a new way to discover short-period 
non-eclipsing binaries with low-mass companions by using highly precise photometri c 



light curves obtained by space missio ns, like CoRoT and Kepler ( IRouan et al. 



Baglin et al. 


2006 


by 


Loeb & Gaudi 



1998 



Borucki et al. 



fl2003f ) and 



2010h. The algorithm, BEER, 



Dased on an idea suggested 



Zucker. Mazeh fc Alexander! (120071 ). searches for the beaming 



effect, sometimes called Doppler boosting, induced by stellar radial motion. This effect 



causes an increase (decrease) of the brig 



from) the observer (IRybicki fc Lightman 



itness of any light source approaching (receding 



19791 ). on the order of 4v r /c, where v r is the radial 



velocity of the source, and c is the velocity of light. Therefore, periodic modulation of the 
stellar velocity due to a companion in a binary orbit will produce a corresponding periodic 
beaming modulation of the stellar photometry. 

For short-period binaries the beaming effect is extremely small, on the order of 100-300 
ppm (parts per million). Therefore the effect has become relevant only recently, when 
CoRoT and Kepler — the two presently operating satellites that search for transiting 



exoplanets, sta rted producing 



high precision ( lAuvergne et al. 



mndreds of thousands of uninterrupted light curves with 



2009 



Koch et al. 



20101 ). 



As predicted, several studies detected the beaming effect in eclipsing binaries and 
transiting pl anets, for which the orbital period was well established from the space-obtained 



light curves ( 



2010 



van Kerkwiik et al 



Bloemen et al. 



2011 



2010 



Rowc et al 



Kipping fc Spiegel 



2010: 



Carter et al 



2010 



Mazeh fc Faigler 



20111 ). Yet, space mission data can yield much 



more. Evidence of the binarity of a st ellar system can be f o und from detecting the beam ing 



effect without any eclipse or transit (ILoeb fc Gaudi 



20031 : 



Zucker. Mazeh fc Alexander 



20071 ). However, the beaming modulation by itself might not be enough to render a star a 
good binary candidate, as periodic modulations could be produced by other effects, stellar 
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variability in particular (lAigrain. Favata &: Gilmordl2004j ). 



The BEER detection algorithm ( IFaigler fc Mazeh 



201 ll ) . therefore, searches for stars 



that show in their light curves a combination of the BEaming effect with two other effects 
induced by the presumed companion — the Ellipsoidal and the Reflection modulation. 



The ellipsoidal variation ((Morris 



19851 ) is due to the t idal d istortion of each component by 



the gravity of its companion (see a review by iMazehl (120081 )). while the reflection/heating 
variation (referred to herein as the reflection modulation) is induced by the luminosity 



of each component that f alls only on the close side of its companion ( Vaz 



1990 



Maxted et al. 



2002 



Harrison et al 



2003 



For et al 



2010 



Reed et al 



1985 



Wilson 



2010|). Detecting 



the beaming effect together with the ellipsoidal and reflection modulations, with the 
expected relative amplitudes and phases in particular, can suggest the presence of a small 
non-transiting companion. 

Just as transit searches, the candidates found by the BEER algorithm have to be 
followed by radial-velocity (RV) observations, in order to confirm the existence of the 
low-mass companion, and to reject other possible interpretations of the photometric 
modulation. 

This paper presents the discovery of the first seven new binaries with low-mass 
secondaries, in the range of 0.07-0.4 M , detected by using the BEER algorithm, and 
confirmed by RV spectroscopic follow-up measurements. Section 2 presents the photometric 
analysis of the Kepler light curves, Section 3 provides the details and results of the RV 
observations, Section 4 summarizes and compares the results of the photometric analysis 
and the RV measurements, and Section 5 discusses the implications of, and conclusions 
from, the findings of this paper. 
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2. Photometric analysis 



We used the publicly available Kepler raw light curves of the QO, Ql and Q2 quarters, 
spanning 138 days. To avoid systematic variations, we ignored all data points within 1 day 
after the beginning of Q2, and all data points within 1 day before, to 3 days after, each of 
the two safe mode events in Q2. We also corrected two systematic jumps at Kepler time 
(JD - 2454833) of 200.32 and 246.19 days. We then appli ed the BEER algo rithm to 14,685 



stars brighter than 13th mag, with Kepler Inpu t Catalog (IBrown et a 



than 3R Q , calculating the BEER periodogram (IFaigler fc Mazeh 



20111 ) radius smaller 



20 111 ) with period range of 



0.5-20 days for each star. Next, we identified the periodograms whose highest peak was at 
least 3 time higher than the next highest one. For these stars we used the peak period to 
estimate the system secondary mass and radius, assuming the periodicity is induced by a 
secondary star. We then selected 25 candidates with secondary mass smaller then O.5M 
and implied albedo smaller then 0.4, suggesting a significant probability for a low-mass 
companion. These candidates were then followed by RV observations, which we describe in 
detail in the next section. In a forthcoming paper we will report on the false alarm cases, 
and analyze the false alarm frequency of our candidates. Here we report on the first seven 
clear detections. 

Table 1 lists for each of th e seven stars its co ordinates, the stellar properties estimates 



from the Kepler Input Catalog (jBrown et al. 



201 ll ) . the photometric periods and amplitudes 



of the three effects found by the BEER algorithm, and the r.m.s. of the data before and 
after subtraction of the BEER model. 



We order the stars according to the detected RV amplitud e 



section. Figure 1 presents the 'cleaned' (IMazeh fc Faigler 



2010 



, presented in the next 



Faigler fc Mazeh 



201 lh 



photometric data of the seven detections, Figure 2 presents the BEER periodograms for the 
detections, and Figure 3 shows the light curves folded with the detected period. In fact, the 
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quality of the Kepler data is so high that the periodic modulation can be seen directly from 
the cleaned data, plotted in Figure 1, even without consulting the BEER periodogram. 

It is interesting to compare the shape of the BEER modulation of the seven candidates, 
presented in Fig 3. In six of them the two peaks, at phase of 0.25 and 0.75, are similar, 



201lh . 



although the latter is somewhat smaller, due to the beaming effect ( jFaigler fc Mazehl 
In one case, K08016222, the second peak completely disappeared, because in this case the 
beaming amplitude is more than three times higher than that of the ellipsoidal, while for 
the rest of the candidates, the ellipsoidal amplitude is significantly higher then the beaming 
amplitude. This is a clear result of the long orbital period and small stellar radius of this 



system, relative to the other syst ems, since the el 



ratio is proportional to ft/I™ kider fe Mazeh 



ipsoidal amplitude to beaming amplitud e 



2011 



Zucker. Mazeh fc Alexander 



20070 . 
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Table 1. Coordinates, magnitudes, stellar properties, and photometric analysis results of 

the seven candidates 





K10848064 


K08016222 


K09512641 


K07254760 


K05263749 


K04577324 


K06370196 


RA 


19:01:21.24 


19:06:48.03 


18:58:39.91 


18:42:28.78 


19:12:59.00 


19:42:35.91 


19:35:00.36 


DEC 


48:16:32.90 


43:48:32.90 


46:08:52.80 


42:49:31.90 


40:26:42.30 


39:38:00.80 


41:47:59.60 


K p a [mag] 


12.13 


11.65 


11.66 


12.04 


11.53 


11.98 


11.97 


R a [Rq] 


1.5 


1.3 


1.7 


1.5 


1.9 


1.3 


2.1 


M b [M Q ] 


1.2 


1.1 


1.2 


1.2 


1.3 


1.2 


1.3 


Photometry results: 
















Period [days] 


3.49 ± 0.01 


5.60 ±0.02 


4.65 ±0.02 


2.66 ±0.01 


3.73 ±0.01 


2.33 ±0.01 


4.23 ±0.01 


Ellipsoidal [ppm] 


201 ± 3 


30 ±2 


172 ± 10 


845 ± 7 


1222 ± 5 


1489 ± 4 


1210 ± 10 


Beaming [ppm] 


118 ±3 


97 ±2 


185 ±5 


356 ±6 


358 ± 5 


436 ±4 


382 ± 7 


Reflection [ppm] 


0±3 


6±2 


36 ±4 


150 ±6 


158 ± 5 


245 ±4 


174 ±7 


Cleaned data r.m.s. [ppm] 


204 


106 


227 


807 


1184 


1408 


1141 


Residuals r.m.s. [ppm] 


128 


68 


113 


268 


200 


168 


328 



a from Kepler Input Catalog 

6 calculated from Kepler Input Catalog log g and R 
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Fig. 1.— The light curves of the seven detections, after outlier removal and long-term 
detrending. Top to bottom: K10848064, K08016222 , K09512641, K07254760, K05263749, 
K04577324, K06370196. For clarity, each light curve was shifted by 5000 ppm relative to 
the previous one. The periodic modulation can be seen in all seven light curves. The light 
curves show several discontinuities: end of QO at day 131, end of Ql around day 166, Q2 
first safe mode event at day 183, and Q2 second safe mode event at day 232. In addition 
there is a single discontinuity at day 155 of the K05263749 light curve. 
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Fig. 2. — The BEER periodograms of the seven candidates. The peak frequency corresponds 
to the suspected orbital period. The periodograms were calculated for the period range of 
0.5-20 days. For clarity, only the period range of 2-20 days is plotted, since no significant 
peak was found for periods smaller than 2 days in any of the periodograms. The periodograms 
are normalized so that the r.m.s. of the 100 noise points on two sides of the peak (50 on 
each side) is set to one. 
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Fig. 3. — Folded cleaned light curves, binned into 100 bins, of the seven candidates. Phase 
zero is when the presumed companion is closest to the observer, while phase 0.5 is when the 
primary is closest to the observer. The errors of each bin represent la estimate of the bin 
average value. The line presents the BEER model. The model residuals are plotted at the 
bottom of each figure. 
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3. Radial velocity observations 

The RV observations were performed between 25 Septe mber 2010 and 15 June 2011 



Fiireszl (120081 )) mounted on 



with the Tillinghast Reflector Echelle Spectrograph (TRES; 
the 1.5-m Tillinghast Reflector at the Fred Lawrence Whipple Observatory operated by the 
Smithsonian Astrophysical Observatory (SAO) on Mount Hopkins in Southern Arizona, 
using the medium fiber at a spectral resolution of 44,000, covering a spectral range from 385 
to 910 nm. Exposures of a Thorium- Argon hollow-cathode lamp immediately before and 
after each exposure were used for wavelength calibration. The spectra were extracted and 
rectified to intensity v s. wavelength using standard procedures developed by Lars Buchhave 



(IBuchhave et al 



2010h 



To derive precise relative radial velocities, we performed a cross-correlation between 
each observed spectrum and a template spectrum constructed by shifting and co-adding all 
the observed spectra. In addition to the template constructed by shifting and co-adding all 
the observed spectra, we also tried using the strongest individual exposure of each object 
as the observed template. The two approaches gave essentially indistinguishable results, 
with slightly better residuals from the orbital fits for the shifted and co-added template. 
We also derived absolute velocities using the library of synthetic templates and found the 
same orbits, although with somewhat larger residuals. 

We did not include spectral orders that were significantly contaminated by telluric 
lines from the Earth's atmosphere, nor did we include the bluest orders with the lowest 
signal-to-noise ratio and a few red orders with known problems. The error of each relative 
velocity was estimated using the standard deviation of the velocities from the 21 individual 
orders, but the velocities themselves were derived by first co-adding the correlation functions 
from the 21 orders, to get a natural weighting of the contribution from each order. 

Using the shifted and co-added template can distort the cross-correlation peak because 
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the noise in each spectrum correlates with the same noise that is still present in the 
averaged template, and therefore can lead to underestimated uncertainties of the velocities. 
To correct this effect we later inflated the uncertainties of the orbital elements (see Xred 
discussion below). 

We used a library of synthetic spectra, calculated by John Lair d for a grid o 



Kurucz mode 



Latham et al 



atmo spheres, using a line list developed by Jon Morse (jCarney et al 



1987 



20021 ). to estimate values for the effective temperature, surface gravity, 
metallicity, and rotational velocity of the seven primaries. This was done by cross-correlating 
each coadded observed template spectrum against a grid of synthetic templates surrounding 
the one that gave the best correlation. Our library of synthetic spectra has a spacing 
of 250 K in effective temperature, T eff ; 0.5 in log surface gravity, logg; 0.5 in the log of 
the metallicity compared to the sun, [m/H]; and has a progressive spacing for rotational 
velocity, v TOt . Because the grid is coarse, we used the correlation peak heights to interpolate 
between grid points, arriving at a more precise classification. Three TRES spectral orders 
overlap with the synthetic spectra, so we performed this cross-correlation and interpolation 
in each order. The mean values, weighted by the cross-correlation peak height in each order, 
and RMS errors are reported in Table 2. Note that because of the degeneracies between T e s, 
logg, and [m/H] in the stellar spectra, correlated systematic errors may dominate. For this 
reason, and based on our experience in other surveys we have inflated the errors by adding 
100 K in T e ff and 0.1 dex in log g and [m/H] in quadrature to the formal order-to-order 
RMS errors. 

The relative velocities were adjusted by a constant offset to a system of absolute 
velocities using observations of the nearby IAU Radial Velocity Standard Star HD 182488, 
whose absolute velocity was assumed to be —21.508 km s" 1 . This adjustment utilized our 
library of synthetic templates, from which we picked the synthetic template that gave 
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the best match to the observations of each star in the spectral order centered on the Mg 
b feature near 518 nm. This approach should avoid the problem of possible template 
mis-match between the various target stars and HD 182488. The uncertainty in the zero 
point of our absolute velocities is probably limited by the uncertainty in the absolute 
velocity of HD 182488, which could be as large as 100 m s _1 . Table 3 lists the radial- velocity 
measurements and their uncertainties. 

For all seven candidates discussed here the first RV measurements showed clear 
variability. We therefore obtained enough RVs to allow orbital solutions completely 
independent of the BEER analysis. To determine the orbital elements of each target, 
independent of the BEER results, we ran a Markov chain Monte Carlo (MCMC) analysis of 
the radial velocities. We adopted values for the epoch (T), period (P), systemic velocity (7), 
orbital semi- amplitude (K), eccentricity (e), and argument of periapse (w) corresponding 
to the median values of the posterior distributions. The errors listed in the tables are those 
corresponding to the lQ th and 84 th percentiles of the posterior distributions. The reported 
error on 7, however, includes contributions both from the formal error from the MCMC 
posterior and from the uncertainty in the TRES absolute zero point offset. 

When the orbit is circular, the epoch reported is T max , the time of maximum velocity, 
and when the orbit is eccentric, we report T per i, the time of periastron passage. In six of 



Table 2. Spectra derived stellar properties of the seven binaries 





K10848064 


K08016222 


K09512641 


K07254760 


K05263749 


K04577324 


K06370196 


T cB [K] 
logp [dex] 
[m/H] [dcx] 
usini [km s _1 ] 


6209 ± 131 
3.68 ±0.16 
-0.24 ±0.11 
14.26 ±0.45 


5919 ± 128 
4.29 ±0.14 
-0.14 ±0.11 
3.37 ±0.24 


6348 ±214 
4.04 ±0.23 
-0.37 ±0.15 
9.47 ±0.24 


6377 ± 133 
4.04 ±0.23 
-0.04 ±0.11 
14.96 ± 0.52 


6328 ± 119 
3.54 ±0.14 
-0.21 ±0.10 
21.50 ±0.16 


6515 ± 145 
3.71 ±0.14 
-0.08 ±0.11 
28.90 ±0.36 


6213 ± 168 
3.91 ±0.13 
-0.26 ±0.16 
20.93 ± 0.45 
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Table 3: Radial velocities of the seven binaries 



Time [HJD-2455000] 


RV [m s -1 ] 


a [m s- 1 } 


Time [HJD-2455000] 


RV [m s _1 ] 


a [m s 1 ] 


K10848064: 


K07254760: 


464.710256 


-8835 


119 


694.814729 


37083 


39 


469.624378 


-17974 


76 


695.834208 


-11894 


46 


488.634305 


-16743 


77 


697.807122 


15727 


61 


489.574441 


-6552 


123 


699.779976 


46148 


83 


490.617670 


-17742 


83 


702.796583 


36412 


47 


498.599389 


-23508 


57 


703.815335 


-11560 


54 


513.639601 


-9108 


170 


704.799619 


39179 


51 


692.888788 


-13130 


124 


705.802918 


14063 


89 


722.832949 


-14423 


43 


K05263749: 


723.863233 


-7328 


56 


694.829973 


20947 


69 


724.857438 


-21205 


28 


696.825053 


13787 


141 


K08016222: 


697.826976 


45095 


68 


465.787007 


-31819 


63 


699.974074 


-12159 


61 


466.613973 


-37012 


39 


701.980887 


34841 


87 


467.723417 


-33482 


96 


705.958593 


23622 


73 


469.613641 


-18312 


42 


722.851562 


10492 


63 


485.603398 


-23367 


37 


723.875349 


45267 


46 


490.579108 


-29588 


63 


724.846275 


10122 


75 


498.591317 


-22720 


47 


K045 77324: 


722.843084 


-21815 


86 


695.982214 


-22723 


166 


K09512641: 


697.845760 


8912 


99 


658.935318 


6324 


22 


703.971684 


38803 


87 


669.855208 


33915 


29 


722.871612 


46652 


99 


693.803705 


34608 


67 


723.895443 


-22140 


103 


694.822698 


16733 


54 


724.869757 


34694 


124 


696.815346 


17503 


50 


726.867639 


5528 


96 


697.816773 


34761 


48 


727.847462 


33003 


91 


698.934501 


27404 


39 


K06370196: 


701.814582 


24838 


40 


694.852291 


20398 


58 


702.805990 


36143 


39 


697.835776 


-13709 


144 


703.826109 


22562 


59 


702.968624 


22603 


41 


704.813460 


6060 


65 


722.861284 


-35084 


114 


705.813161 


11855 


43 


723.884266 
724.826295 
726.857187 
727.797595 


17354 
8740 
-47080 
4696 


76 
72 
95 
132 



seven cases, either the orbital phase coverage was not sufficient to adequately constrain the 
eccentricity or e was statistically indistinguishable from zero. In these cases, we fixed e = 
and reran the MCMC chains, adopting T max , P, 7, and K from this solution. In one case, 
K080 16222, the orbital phase coverage is good and e is significantly non-zero. 

Figure 4 shows the RV follow-up measurements for each of the seven binaries, folded 
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Fig. 4. — The RV measurements of the seven stars, folded at the derived orbital period. 
The line presents the orbital RV model. The residuals are plotted at the bottom of each 
figure. Note the different scale of the upper and lower panel of each star. The error bars are 
too small to be seen in the upper panels. 
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with the period found, and Table 4 lists the derived orbital elements. The table also lists 
Xred' the reduced x 2 °f the model, and the time span of the observations. For two binaries 
the derived x 2 e d value is close to unity, as expected, but for the others its value is relatively 
large. This could indicate either that for those binaries our radial-velocity uncertainties 
are underestimated, or that our radial-velocity model is too simple, due to some stellar 
noise, for example. In order to get more realistic uncertainties for the model elements, we 
inflated the parameter uncertainties of each target by its \/xl e di which is equivalent to 
inflating the RV errors of that star by the same factor. The resulting orbital model elements 
uncertainties are listed in Table 4. 



4. Results 

Table 5 lists for each of the seven newly discovered binaries the period derived 
from the photometry, the calculated abeam, and the expected RV semi-amplitude, Kb eam , 
derived from abeam and the photometric beaming amplitude. The abeam factor includes 
one component that originates from the fact that the stellar spectrum is Doppler shifted 
relative to the observed band. To estimate this factor fo r each of the seven detec ted 
binaries we numerically shifted spectra from the library of ICastelli &: Kuruca ( 12004 ) models 



that were close to the estimated temperature, metallicity and gravity of each of the seven 
stars. The values adopted were derived by interpolation of the abeam values between the 
available models of the library. The abeam uncertainties were estimated by calculating the 
interpolated abeam values within the T eS , log g and [m/H] error ranges. The error on the 
expected Kb ea m was estimated by combining the photometric beaming amplitude error and 
the abeam error. The table then reports the number of RV measurements, their derived 
RV period and semi-amplitude, and the minimal secondary mass, up to sinz. For all cases 
we independently derived the period of the RV modulation, and found it to be consistent 
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Table 4. Orbital model elements of the seven binaries 





K10848064 


K08016222 


K09512641 


K07254760 


K05263749 


K04577324 


K06370196 


T max - 2455000 [HJD] 


465.1289 


464.6288° 


642.3690 


702.4395 


466.7071 


466.7005 


698.8374 




±0.0060 


±0.0029 


±0.0052 


±0.0126 


±0.0124 


±0.0079 


±0.0035 


P [days] 


3.49318 


5.60864 


4.64588 


2.65642 


3.72665 


2.328663 


4.23371 




±0.00099 


±0.00017 


±0.00044 


±0.00068 


±0.00019 


±0.000070 


±0.00067 


7 [km s _1 ] 


-15.670 


-28.078 


20.518 


17.092 


13.862 


11.325 


-17.168 




±0.219 


±0.048 


±0.110 


±0.177 


±0.060 


±0.119 


±0.454 


K [km s- 1 ] 


9.107 


9.495 


15.519 


29.024 


31.428 


35.316 


40.222 




±0.073 


±0.018 


±0.023 


±0.061 


±0.040 


±0.043 


±0.131 


e 


(fixed) 


0.0439 ± 0.0022 


(fixed) 


(fixed) 


(fixed) 


(fixed) 


(fixed) 


oj [deg] 




36.2 ± 2.6 












X 2 A 
"■red 


11.0 


1.0 


4.3 


11.6 


2.6 


1.2 


12.6 


span [days] 


260.1 


257.1 


46.9 


11.0 


30.0 


31.9 


32.9 



a for K08016222 the T max value is the time of periastron passage 
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with the photometric period, indicating that the orbital period was reliably derived by the 
BEER algorithm, solely from the photometric data. 

In six of the binaries the eccentricity was too small to be derived significantly, so 
we assumed circular orbits. Because these are short-period stellar binaries, the expected 
circularization timescale is short, so finding in most cases that e = is consistent with our 
expectations. For K08016222 we find e = 0.0439 ± 0.0022. Interestingly, this is the binary 
with the longest period, so its lifetime might have been too short to achieve circularization 
( IMathieu fc Mazehlll988[ ). 

Out of the seven binaries, the measured RV amplitudes of five cases were consistent 
with those predicted by the photometric analysis. For the other two stars, K080 16222 and 
K06370196, the predicted amplitudes were 24% smaller than the observed ones. This could 
be due to underestimation of the photometric amplitude. Another possible explanation may 
be an inaccurate translation of the photometric amplitude to the expected RV amplitude, 
which depends on the assumed stellar spectral type. We need more confirmed binaries to 
understand this effect. 



5. Discussion 

The RV observations presented here demonstrate the ability of the BEER algorithm to 
discover short-period binaries with minimum secondaries mass in the range of 0.07-0.4 M & 
in the publicly available Kepler data. 

The original goal of the Kepler and CoRoT missions was to search for transiting 
planets. Such projects are limited to planets with orbital inclinations cl ose to 90°. The 



serendipitous discoveries of eclipsing binaries in the Kepler photometry (jPrsa et al 



2011) 



are suffering from the same limitation. The BEER algorithm, on the other hand, is 



- 19 - 



Table 5. Derived photometric RV period and semi-amplitude together with RV 
observations period and semi-amplitude for each of the seven binaries 





K10848064 


K08016222 


K09512641 


K07254760 


K05263749 


K04577324 


K06370196 


Photometry results: 
Period [days] 

^beam 

K bearn [km s- 1 ] 


3.49 ±0.01 
0.944 ± 0.025 
9.37 ±0.34 


5.60 ±0.02 
1.012 ±0.023 
7.19 ± 0.22 


4.65 ±0.02 
0.912 ± 0.036 
15.21 ±0.72 


2.66 ±0.01 
0.921 ±0.024 
28.97 ±0.89 


3.73 ± 0.01 
0.921 ±0.022 
29.14 ±0.81 


2.33 ±0.01 
0.887 ±0.025 
36.86 ± 1.10 


4.23 ±0.01 
0.936 ± 0.030 
30.61 ± 1.12 


RV results: 
N obs 


11 


8 


12 


8 


9 


8 


8 


Period [days] 


3.49318 
±0.00099 


5.60864 
±0.00017 


4.64588 
±0.00044 


2.65642 
±0.00068 


3.72665 
±0.00019 


2.328663 
±0.000070 


4.23371 
±0.00067 


K RV [km s _1 ] 


9.107 
±0.073 


9.495 
±0.018 


15.519 
±0.023 


29.024 
±0.061 


31.428 
±0.040 


35.316 
±0.043 


40.223 
±0.131 


Minimum secondary 


76 ±5 


90 ±6 


147 ± 10 


222 ± 15 


279 ± 19 


253 ± 17 


376 ± 25 



mass [M Jnp ] 
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searching for non-transiting companions, and therefore can detect many more systems with 
much lower inclination angles. Searching with BEER is effectively equivalent to performing 
an RV survey that is not limited to nearly face-on inclinations. Applying the BEER 
algorithm to the hundreds of thousands of already available light curves of Kepler and 
CoRoT is like performing an RV survey of a huge sample that is composed of these stars. 

Therefore, we expect BEER to discover many hundreds of new binaries with short 
periods. Furthermore, whereas in RV studies the actual mass of the companion depends on 
the unknown inclination angle, detecting both the ellipsoidal and the beaming effects will 
enable BEER to derive, or at least estima t e, the mass of the small companion in certain 



cases. As pointed out by 



Faigler fc Mazehl ( 120111 ). this can become possible because the two 



effects have different dependencies on the orbital inclination, and therefore the derived ratio 
of the amplitudes of the two effects can, in principle, remove the degeneracy between the 
secondary mass and the inclination. 

Obviously, at this stage of the BEER search, detecting a candidate is not enough - 
the candidates have to be confirmed by follow-up RV observations. However, when we 
accumulate enough observations we will be able to estimate the false alarm probability, 
which might be a function of the amplitude of the photometric modulation and the stellar 
mass, radius and temperature. Therefore, we will be able to derive the statistical features 
of the short-period binaries without confirming each detection with RV observations. 



Faigler fc Mazeh 



he seven cases presented here were based on the Kepler Q0-Q2 data. 
( 1201 ll ) suggested that once the full Kepler dataset is available, we should be able to detect 
brown-dwarf secondaries and even massive planets. Moreover, the other stellar modulations 
that contribute now to the false alarm frequency are not expected to be so stable on time 
scales of years, whereas the three BEER effects are strictly periodic and stable. Therefore, 
we expect the false alarm frequency to decrease when we have access to longer data sets. 
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The unprecedentedly large sample size and data quality, together with a knowledge of the 
false alarm probability, could serve as a tool to study accurately the frequency of low-mass 
second aries in short-period b i naries on the high- and low-mass ends of the brown-dwarf 



desert (IRaghavan et al. 



2010 



Udry 



2010 



Sahlmann et al. 



2010). 
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